The lowest scattering state of one-dimensional Bose gas with attractive interactions 
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We investigate the lowest scattering state of one-dimensional Bose gas with attractive interactions 
trapped in a hard wall trap. By solving the Bethe ansatz equation numerically we determine the 
full energy spectrum and the exact wave function for different attractive interaction parameters. 
The resultant density distribution, momentum distribution, reduced one body density matrix and 
two body correlation show that the decreased attractive interaction induces rich density profiles and 
specific correlation properties in the weakly attractive Bose gas. 
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I. INTRODUCTION 

With the rapid experimental progress the ultracold 
atomic gases have offered a popular platform to investi- 
gate the strongly correlated one dimensional many-body 
systems [T-^ for their high controllability and tunabil- 
ity. One dimensional quantum gases can be realized with 
strong anisotropic magnetic trap or two dimensional op- 
tical lattice and be described theoretically by an 
effective one dimensional (fD) model Jhf Q]. In addi- 
tion the effective ID interaction can be tuned from the 
strongly attractive to the strongly repulsive interacting 
regime via the magnetic Feshbach resonance or confine- 
ment induced resonance. Not only the strongly inter- 
acting Tonks- Girardeau (TG) gases [ll| but also two 
counter-intuitive examples, i.e., the repulsively bound 
atom pairs fl^l and the super Tonks- Girardeau (STG) 
gas-like phase of the attractive atomic gas [T3l - fl5| , have 
been realized by a sudden quench of interaction from the 
strong repulsion to the strong attraction or vice versa, 
both of which are hard to realize in the traditional con- 
densed matter physics and have no analog in solid state 
systems. It has been displayed that they are stable ex- 
cited states and the stability could be understood from 
the quench dynamics of the ID integrable quantum gas 
p^ . So far, the STG gas has attracted intensive theo- 
retical studies from various aspects [T7l - [2lj . 

The experimental realization of STG gases and repul- 
sively bound atom pairs open the door to stable highly 
excited quantum many-body phases. This also offers us 
a method to search for exotic quantum phases in ID 
many-body systems. It has been predicted that such 
stable excited states can be prepared in optical lattice 
via sudden quantum quench |17| and the effective su- 
per Tonks- Girardeau gases can be realized via strongly 
attractive one-dimensional Fermi gases [l^- Theoreti- 
cally the strongly interacting Bose gases cannot be well 



described in the mean-field theory due to the large quan- 
tum fluctuation in ID system and one has to resort to 
non-perturbation methods. Many methods such as Bethe 
ansatz 22 24|, Bosonization method [25], exact diag- 
onalization [26|, l27| , Bose- Fermi mapping method (BFM) 
[28j and multi-configuration Hartree theory [29| were 
used to investigate the ID quantum gases. It was shown 
that with the increase in repulsion strength the ground 
state density distribution of ID Bose gases continuously 
evolves from a Gaussian-like distribution to a multi-peak 
structure while the momentum distribution remains the 
single peak structure of bosonic atoms [l^, [H, [13, [13] ■ 

By a sudden quench of interaction from strong repul- 
sion to stron g at traction a TG gas shall transfer into a 
STG gas Tsl Il6l. [Toj . Because the quantum gas is very 
weakly coupled with the environment its energy dissi- 
pation is ignorable and this highly excited state shall 
be stable. By decreasing the attractive interaction of 
STG gas we can investigate the properties of the lowest 
scattering state for the Bose gas with the change of at- 
tractive interaction from strongly to weakly interacting 
regime. In this work, by numerically solving the Bethe 
ansatz equation we obtain exact wave function of the 
lowest scattering state of ID Bose gas trapped in a hard 
wall potential in the above interacting regime. We will 
focus on the density distribution, reduced one body den- 
sity matrix (ROBDM) and two body correlation in the 
full attractive interaction regime. 

The present paper is organized as follows. Section II 
is devoted to the description of our model and Bethe 
Ansatz method. Section III will give the density profiles, 
ROBDM and two body correlation in the full attractive 
interacting regime. A summary is given in the last sec- 
tion. 



II. MODEL AND METHOD 
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We consider N interacting Bose atoms of mass m con- 
fined in a hard wall box of length L, which is described 
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by the Hamiltonian 
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H = - 



N ^2 

2—1 ^ 



l<i<j<Ar 



(1) 



Here the natural unit is used h — 2m = 1 and c ~ 
TTigio/^^ is an interaction constant dependent on the ef- 
fective ID interaction strength giD, which can be tuned 
continuously from the strong attraction to the strong re- 
pulsion by Feshbach resonance or confinement induced 
resonance. This model can be solved exactly by the Bethe 
ansatz method p}]| . The many-particle wave function 
shall be formulated as the following general form 

Q 

XipQ {Xg,,Xq^,- ■ ■ ,Xg^) (2) 

with 
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where Q = (gi,g2,-- - ,9jv) and P (pi,P2,-- - ,Vn) 
are one of the permutations of I,-- - ,-/V, respec- 
tively, A{Q,rP) is the abbreviation of the coefficient 
^i(li,(l2,- ■ ■ ,qN;rp^Pi,rp^P2,- ■ ■ ,rp^PN) to be deter- 
mined self-consistently, and the summation J^p (Sq) 
is done for all of them. Here rj — ± indicate that the 
particles move toward the right or the left, 9{x — y) is 
the step function and the parameters {kj} are known as 
quasi-momenta. In the following evaluation the length L 
will be taken to be unity unless otherwise specified. For 
Bosons the wave function should follow the symmetry of 
exchange, so the present problem is simplified into the 
solution of 

H(pi,...,N (xi,--- , xn) = Eipi,,,,,N [xi,--- , xn) (4) 

< xn < L with the 



in the region of < xi < 2:2 < 
open boundary condition 



0, 



.N I 



0. 



For simplicity we shall ignore the subscript in 
V^i,...,N {xi, ■ ■ ■ ,xn)- The wave function in other region 
Q can be obtained by the exchange symmetry of Bose 
wave function. 

With some algebraic calculation, the wave function has 
the following explicit form 



(p{xi,X2, ■ ■ ■ ,Xn) 
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FIG. 1: (color online) Energy levels of Bose gas with N — 2 
for different interaction c. 
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Here Ep= ±1 denote sign factors associated with 
even(odd) permutations of P. The quasi-momenta are 
determined by numerically solving the Bethe ansatz 
equations and the total wave function is given by Eq.Q 
through ip{xi, • • • , Xn) under the restriction of exchange 
symmetry. 

Under the open boundary condition, we can obtain 
Bethe-ansatz equations 



J2kiL 
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kj + ki + ic kj — ki + ic 



ki + ki — ic k-i — ki — ic 
i=i(i^j) -J ' ^ ' 

whose logarithmic forms are formulated as 



kjL = TLji: 



N 

E 



kj ~\~ ki kj ~ ki 
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X Yl ^'^^[Kj^j - 2^^P,Pj 

l<j<N \ l<] 



exp (ikpj^L) sin {kp-^xi) 



sin (fcp„ {L - Xn)) 



(5) 

Here rij is a set of integers to determine the eigenstates 
and for the ground state rij — j {1 < j < N). The energy 

of the system is E = J2f=i ^j- 

For the repulsive interaction (c > 0) the solutions of 
Eq. ([5]) kj are real, while for the attractive interac- 
tion (c < 0) its solutions might be either real or com- 
plex. The real solutions correspond to scattering states, 
which can be obtained by solving Eq. We de- 

termine the quantum number Uj in the limit of strong 
interaction (c — >■ ±00) or the limit of zero interaction 
(c ^ 0). In the former situation the exact wave function 
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can be constructed by means of the Bose-Fermi map- 
ping method from the wave function of N free Fermions 
[19l |. i.e., the eigenfunction of single particle in the hard 
wall trap sin^jirx/L) with j being integers. For the 
ground state j — 1, - ■ ■ ,N and the excited states are ob- 
tained when some j are replaced by the integers greater 
than N. In the ground state of noninteracting limit all 
Bose atoms shall condense into the ground state of sin- 
gle particle, which gives the exact many body function 
as (l){xi, • ■ • ,a;jv) = Yij siniirxj / L). The wave function 
of excited states in noninteracting limit, however, takes 
the form of • • • ,xn) = PY\j sin^rijirxj / L), where 
P is an operator preserving the exchange symmetry of 
the many body wave function. By comparing these ex- 
act wave function with those obtained from Eq. ([S]), the 
quantum number rij can be obtained. It is convenient to 
use the strong interacting limit here. As c approximates 
to ±00 all terms in the summation of Eq. ([5]) vanish 
such that the quasimomenta kj = Ujir/L. Comparing 
the above exact solution we have nj — j {j — 1, - ■ ■ ,N) 
for the ground state and the excited states correspond to 
Uj > j. After deciding the quantum number it is easy 
to obtain the solution in the limit of c — ±0. As c ap- 
proximate 0^ we have kj = (uj — j + 1)ti/ L and as c 
approximate 0^ we have kj — {rij + j — 1)tt/L. 

The complex solutions correspond to bound states, 
which can be obtained by assuming the solutions of com- 
plex form kj = Uj +iAj and solving the set of equations of 
Uj and Aj. For example, the complex solutions of Bethe 
ansatz equations for the system of = 2 are assumed as 
kj — k/2 + iA and the Bethe ansatz equations take the 
formulation of 
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exp[2AL] = (-1)" 
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where the integer n > 2 is quantum number and for the 
ground state n = 2. As an example, we display the full 
energy spectrum for A^ = 2 in Fig. 1. The scattering 
states are denoted by solid lines and the bound states 
are denoted by dashed lines. It is shown that the ground 
state for repulsive case is scattering state and that for 
attractive case is bound state. The excited state of at- 
tractive Bose gas in the scattering state has been realized 
experimentally in Ref. [Tsj . By tuning the ID interac- 
tion constant the Bose atoms evolve from the weakly in- 
teracting Thomas-Fermi regime to the strongly repulsive 
TG regime in which the TG gas was realized, and then 
by quenching the interaction from the strong repulsion 
to the strong attraction the scattering states of attrac- 
tive Bose gas, i.e., the STG gas, was realized. Starting 
from a stable STG gas, we can further investigate the 
crossover behavior of the STG gas when one decreases 
the attractive interaction very slowly to the very weak 
limit. Through an adiabatically slow change of the at- 
tractive interaction, the lowest scattering state of attrac- 
tive ID Bose gas in the whole attractive regime could be 
reached. 
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FIG. 2: (color online) The quasimomentum distribution and 
the corresponding energy for the gas with A'^ = 6 and c — 
-1000, -100, -30, -10, (from below to above). 



In Fig. 2 we display the quasimomentum distribution 
and the corresponding energy (longitudinal axis) for the 
lowest scattering state of the attractive Bose gases with 
A^ = 6 in the full attractive interacting regime. Accord- 
ing to Eq. dSJ in the limit of STG (c — > —00) the solutions 
of Bethe ansatz equation are equal to those of TG gas 
(c +00), i.e., kj = j-k/L (j = 1, ■ • • , A^). We find that 
with the decrease of attractive interaction strength, kj 
tend to distribute with larger and larger space between 
them although the lowest quasimomentum fci increase 
first and then decrease. In the limit of noninteracting 



limit c — >■ we have kj — (2j — l)7r/L (j = 1, 



,N) 



and the space between two neighbor kj has evolved from 
tt/L (c = —00) to 27r/L (c = 0). The energy of the lowest 
scattering state increases with the decrease of attractive 
interaction. 



III. PROPERTIES OF THE LOWEST 
SCATTERING STATES FOR ATTRACTIVE 
BOSE GASES 

In terms of the lowest scattering state wave function 
^(a;i,--- ^xjq) the important quantity in one dimen- 
sional interacting many-body system, the ROBDM, can 
be formulated as 



p{x, x) 

_ N dx2 ■ ■ ■ dxN^* {x,X2, - ■ ■ , xn) {x' ,X2, ■ • • , xn) 
dxi ---dxN I* (a;i,a;2, • • • ,XN)f 

The diagonal part of ROBDM gives the expectation val- 
ues of density distribution p(x) = p{x,x')\x=x' , and the 
off-diagonal part gives information of coherent properties 
of the gas. The Fourier transformation of p{x, x') gives 
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FIG. 5: (color online) Momentum distribution for the gas 
with iV = 4. 



FIG. 3: The reduced one body density matrix for the gas with 
iV = 4 for c = -1000, -10, -1, (from the left above to the 
right below). 




FIG. 4: (color online) Density distribution for the gas with 
N = 4. 



the momentum distribution 



n{k) 



1 



dx / dxV(x,x')e-*'=("-"'). (6) 







Fig. 3 shows the ROBDM for the lowest scattering 
state of attractive Bose gas with = 4. It deserves 
to notice that for all attractively interacting strengthes 
there exists a strong enhancement of the diagonal con- 
tribution p{x,x') along the line x = x' . In the limit of 
strong attraction, the system displays the same behav- 
ior as TG gases that p{x,x') reduces rapidly as \x — x'\ 
increases. For weaker attraction the off-diagonal part of 
ROBDM shall increase gradually and we have approxi- 
mately p{x, L — x) ~ p{x, x) when the attraction is weak 



enough. The density distributions of the lowest scatter- 
ing state for different attractive interacting strength are 
displayed in Fig. 4. In the strongly attractive limit, it is 
shown that the density profile of the STG gas of A'' atoms 
exhibits the Fermi-like shell structure of A^-peak similar 
to the density profile of the TG gas. This is due to the 
fact that the STG state in the limit of c ^ — oo and the 
TG state in the limit of c — ?> oo are actually identical. 
As the attraction decreases, the density distribution de- 
viates the Fermi-like distribution and the shell structure 
oscillates more and more dramatically. In the limit of c 
approaching 0~ there appear 2N — 1 peaks in the density 
profile. The atoms tend to populate at the center of the 
trap with the most probability and the density distribu- 
tion displays an obvious peak in the center and oscillates 
in the region away from the center. In the limit of the 
strong attraction atoms populate in the N lowest eigen- 
states of single particle so the density profiles show the 
structure of A^-peak. As the decrease of attraction atoms 
populate at higher eigenstates and the peak number shall 
increase. As c — ?> 0~, atoms distribute at the A^ lowest 
odd states of single particle such that the density profile 
exhibits 2N — 1 peaks. 

The momentum distributions for the lowest scattering 
state are shown in Fig. 5. For strong attractive inter- 
action, the atoms accumulate in the central regime close 
to zero momentum and the population distributions de- 
crease rapidly for large momentum, which reflects the 
statistics of bosonic atoms. With the decrease of at- 
traction the atoms distribute more extensively in higher 
momentum region and in the strong but finite attrac- 
tive interaction (e.g. c = —50) the atoms distribute 
widely in the momentum space without an obvious zero- 
momentum peak. For even weaker attractive interaction, 
the probability of atoms locating at the zero momen- 
tum increases and the momentum distribution develops a 
prominent peak at zero point k = 0. Away from the cen- 
tral peak, the momentum distribution starts to display 
shell structure, which however exhibits distinct feature 
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FIG. 6: Two-body correlation for c = -1000, -10, -5, (from 
the left above to the right below). 

from the momentum distribution of Fermi gas. It is be- 
cause the lowest scattering state is an excited state that 
the Bose atoms populate in higher momentum region. 

It is also interesting to study the two body correlation 
function defined as 

92ix,x') 

_ N{N -1) dX3---dXN\'i! {X,X',X3,--- ,XN)f 
Jq dXi---dXN\'b {X1,X2,--- ,XN)f 

which denotes the probability that one measurement will 

find an atom at the point x and the other one at the 
point x'. In Fig. 6 we display the two body correlation 
of lowest scattering state for the Bose gas with N = A. It 
turns out that two atoms with strong attractive interac- 
tion would try to avoid each other and try to keep away 



from each other in certain distance, while as the decrease 
in attraction the probability of finding two atoms in the 
adjacent region increase and arrive at the maximum as 
c-)- 0. 



IV. SUMMARY 

In conclusion we have investigated the lowest scatter- 
ing state of Bose gas in the full attractive interacting 
regime with Bethe ansatz method. By solving the Bethe 
ansatz equations numerically the exact wave functions of 
the lowest scattering state were determined. Based on 
the wave function we obtain the ROBDM, density pro- 
file, momentum distribution and two-body correlation. It 
is shown that in the STG hmit the ROBDM of the low- 
est scattering state exhibit the same behavior as TG gas. 
With the decrease of attractive interaction the density 
distribution evolves from a iV-peak shell structure to a 
(2A/^ — l)-peak one and Bose atoms are located at the cen- 
ter of the trap with the most probability. The momentum 
distribution manifests the nature of Bose statistics in the 
STG limit with an obvious zero-momentum peak. When 
the attraction keeps on deviating the STG limit the mo- 
mentum distribution spreads more widely although the 
most probable position at which the atoms populate is 
still near the region of zero momentum. The change of 
the two body correlation function with the decrease in 
the attractive interaction is also discussed. 
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